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APPA R I EMENT FIN D'IMAGES STEREOSCOPIQUES ET INSTRUMENT 
DEDIE AVEC UN FAIBLE COEFFICIENT STEREOSCQPIQUE 

La presents invention concern e le domaine de la stereoscopie. En 
particulier, ['invention concerne un precede et un ensemble d'acquisition et 
d'appariement fin d'un couple stereoscopique d'images. 

C'est-a-dire que I'invention concerne un procede et un systdme permettant 
I'acquisition et la mise en correspondance des points d'une premiere Image 
avec les points d'une deuxieme Image formant un couple stereoscopique 
d'images. 

PRESENTA TION GENERALE DE L'ART ANTERIEUR 

La ster6oscople est un procede donnant I'impresslon de relief a partir 
d'un couple d'images 2D (ou deux dimensions) representant une scene ayant 
ete acquise avec des points de vue dlff^rents. 

Tel qu'illustre a la figure 3, deux images formant un couple 
stereoscopique d'images sent acquises a I'aide de deux capteurs CCD (Charge 
Couple Device) (matrices ou barrettes) 15, 16 situes dans le meme plan focal 
19' et symetrlques par rapport a la droite A-A' passant par le centre de la scene 
observee et perpendiculaire au plan focal 19'. 

Ces capteurs CCD vont permettre I'acquisition de deux images de la 
scene observee situee a une certaine distance 18 des capteurs CCD. 

Du fait de I'ecartement 19 des deux capteurs CCD qui les ont acquises, 
ces deux images ne sont pas sans relation. Ainsi, une scene verra ses images 
decalees d'une image du couple a I'autre. Ce decalage est appele disparite. 

Le systeme stereoscopique le plus eIabor§ a I'heure actuelle est le 
systeme visuel humain. L'homme dispose de deux yeux situes dans le meme 
plan focal comprenant chacun un cristallin (systeme optique) et une retine 
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(capteur CCD). Ces yeux lui permettent de percevoir son environnement de 
deux points de vue legerement differents. La fusion mentale de ces deux vues 
lui permet d'engendrer une image spatiale (3D) de son environnement. 

Les methodes de traitement damages stereoscopiques essaient, par des 
calculs mathematiques, de modeliser ce que notre cerveau fait sans difTiculte, 
Le principe consiste dans un premier temps a mettre en correspondance les 
points des images du couple stereoscopique d'images. Dans un deuxieme 
temps, on reconstruit les coordonnees 3D des points de la scene observee en 
calculant ['intersection des droites passant par les points apparies et par les 
centres optiques des capteurs CCD. 

La principale difficulte rencontree lors de la mise en place d'une telle 
methode est justement de mettre un a un en correspondance les points des 
images du couple. En effet, les deux images du couple n'etant pas prises sous 
la meme incidence, un point donne de la scene dont ia position dans la premiere 
image est donnee par les cordonnees (X1,Y1), aura pour coordonnees (X2,Y2) 
dans la deuxieme image, avec X1 9^ X2 et Y1 ^ Y2, 

Cet ecart de position entre les projections sur les capteurs CCD d'un 
point donne de la scene observee dont les deux images se ferment sur des 
points geometriquement non correspondants est appele disparite geometrique. 
Dans la suite, le terme « d'ecarts de position » sera prefere au terme 
« disparite » bren qu'ils aient tous deux la meme signification, 

II a deja ete propose un certain nombre de methodes permettant 
Tappariement des points de 2 images prises sous 2 angles differents. Ces 
methodes utilisent souvent le principe de correlation pour faire correspondre les 
points d'une image a Tautre et done determiner les ecarts de position, 

Le principe de la correlation est fonde sur la mesure d'une ressemblance 
locale entre deux images. Cette ressemblance locale s'effectue en introduisant 
des fenetres ponderees centrees sur ies voisinages homologues au sens du 
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posiflonnement g^ometrique dans ['Image. Ces fenetres ponderees sont 
appelees fenetres de correlation. 

Comme presents sur la figure 1, le procede consiste a se donner una 
fenetre 3 de correlation centree sur !e point 40 etudie de la premiere image 1, et 
a rechercher son homologue 41 au sens radiometrique sur la deuxieme Image 
2. Cette operation est realisee par deplacement de la fenefre 3 de correlation 
dans une fenetre plus importante appelee aire de recherche 4. L'aire de 
recherche 4 est centree sur I'homologue geometrique estlme 42 du point 
courant 40 de la premiere image. Le tableau des valeurs de con-elation 
obtenues constitue la nappe de correlation. On retient alors I'ecart de position 
pour iequel la valeur de correlation est maximale. 

Lorsque I'on effectue Tappariement des points d'un couple 
stereoscopique d'images par correlation, la premiere image sert de reference, et 
I'on cherche a retrouver, pour chaque point de la premiere image, son 
homologue radiometrique dans la deuxidme Image. 

Cette operation est realisee successlvement sur un sous ensemble ou 
sur tous les points de la premiere Image, et le point de la premiere image etudie 
a un instant donne sera dans la suite appele point courant. Chaque point de la 
premiere image peut etre considere comme un signal que I'on cherche a 
retrouver dans la deuxieme image par correlation. Un point 40 de la premiere 
image et un point 41 de la deuxieme image sont des homologues 
radlom^triques s'ils correspondent au m§me point de la scene representee dans 
les Images du couple stereoscopique. 

A partir de deux images 2D formant un couple stereoscopique d'images 
il est possible, une fois les points des deux images du couple apparies. de 
detemiiner la coordonnee en 2 d'un point de la scene ayant pour coordonnees 
(X1, Y1) dans la premiere image, et (X2. Y2) dans la deuxieme image du 
couple. 
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En ['absence de perturbations (microvibrations par exemple) la precision 
sur la determination de la troisieme coordonnee d'un point de la scene aussi 
connue sous le nom de precision altimetrique est proche de : 

dz = ds/{blh)^ 

ou -ds est la precision 17 de la restitution de la disparite, et 

- b/h est le coefficient stereoscopique avec 

- b est la distance 19 entre les deux capteurs CCD situes dans le meme 
plan focal, et 

- h la distance 1 8 entre les capteurs CCD et la scene visee, 

Les methodes qui utilisent le principe de correlation pour apparier chaque 
point de la premiere image avec son homoiogue radiometrique dans la 
deuxieme image ne permettent pas d'obtenir une bonne precision sur la 
restitution de la disparite {ds). En effet, le defaut majeur du precede de 
correlation est un phenomene dit d'adherence aux zones les plus singulieres 
des images du couple stereoscopique, comme par exemple les contours. Le 
phenomene d'adherence est cause notamment par la taille de la fenetre de 
correlation utilisee pour appareiller les points du couple damages. 

Ce phenomene d'adherence est mis en evidence a la figure 2 sur un 
decalage 6 sinusoTdal de revolution de periode 35 et d'amplitude +/- 0.5, la 
mesure 7 de ce decalage 6 par correlation avec une fenetre de correlation de 
taille 15x15 de type prolate et la premiere image 9. La mesure 7 de ce decalage 
par correlation fait apparaTtre le defaut majeur de la correlation au niveau des 
contours 9, 10, 11, 12, 13, de Timage 7. 

Par consequent, I'utillsation du principe de correlation pour determiner les 
ecarts de position entre les points de deux images stereoscopiques necessite 
d'avoir un fort coefficient stereoscopique afin de diminuer lineairement Terreur 
altimetrique, ou erreur dans la determination de la troisieme coordonnee Z d'un 
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point de la scene representee dans le couple d'images. Inversement les faibles 
coefficients stereoscopiques augmentent lineairement I'erreur altimetrique. 

Un des buts de la presente invention est de fournir une methode 
d'appariement d'un couple stereoscopique d'images ne presentant pas les 
inconvenients de ia methode precedemment decrite, et ainsi permettre la mise 
en place d'une methode permettant I'appariement precis des points d'un couple 
stereoscopique dimages pour de faibles coefficients stereoscopiques. 

Cet appariement des points d'un couple stereoscopique d'images est 
realise par Telaboration de cartes des ecarts de position, chaque point de ces 
cartes representant l'6cart de position entre un point de la premiere image et 
son homologue dans la deuxieme image. 

Comme ii a ete precedemment mentionne, les methodes d*appariement 
d'un couple stereoscopique d'images basees sur le princlpe de correlation ne 
permettent d'obtenir un appariement precis des points que pour un fort 
coefficient stereoscopique. 

Le coefficient stereoscopique (b/h) est le rapport entre I'ecartement 19 
entre les deux capteurs CCD 15, 16 situes dans le meme plan focal 19' (cf, 
figure 3), et la distance 18 entre les capteurs CCD et la scene visee. 

La distance 18 entre les capteurs CCD et la scene observee depend de 
I'application. Ainsi, suivant rapplication, i'ecartement 19 entre les capteurs CCD 
15, 16 doit etre prevu suffisant pour avoir un fort coefficient stereoscopique. 

Les systemes d'acquisition et d'appariement de couple stereoscopique 
d'images comprennent un systeme d'acquisition et un systeme de traitement. 
Ces deux systemes sont en general distants et communiquent par des moyens 
de communication avec ou sans fil. 

Les systemes de traitement permettent Tappariement d'un couple 
stereoscopique dimages. Ces systemes mettent en oeuvre les methodes 
d'appariement d'un couple stereoscopique. Dans le cas d'un systeme de 
traitement mettant en oeuvre une methode basee sur le principe de correlation, il 
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est done n§cessaire d'avoir un fort coefficient stereoscopique pour un 
appariement precis des images d'un couple. 

Les systemes d'acquisition de couple stereoscopique d'images 
fournissant aux systemes de traltement les images a traiter doivent done etre 
congus de sorte qu'ils respectent cette condition (coefficient stereoscopique fort 
du couple d'images). 

Dans le cas d'un systeme spatial d'acquisition d'un couple 
stereoscopique d'images, la distance 18 entre le systeme d'acquisition et la 
scene observee est tres grande. 

Par consequent, I'ecartement 19 entre les capteurs CCD du systeme 
d'acquisition doit etre tres important pour avoir un fort coefficient 
stereoscopique. L'importance de cet §cartement Impllque que les systemes 
spatiaux d'acquisition de couple stereoscopique d'images comprennent deux 
Instruments optiques 15', 16' (satellite) comprenant chacun un capteur CCD 15, 
16. 

La presence, dans les systemes spatiaux d'acquisition, de deux 
Instruments optiques 15', 16' comprenant chacun un capteur CCD 15, 16 induit 
une plus grande complexite dans I'appariement d'image, et un cout plus 
important pour la conception et la mise sur orblte de tels systemes spatiaux 
d'acquisition. 

En vision de pres (appareils pfiotos par exemple), I'acquisition d'images 
stereoscopiques necessite la presence de deux systemes optiques dans le 
systeme d'acquisition. 

Un autre but de la presente invention est de fournir un ensemble 
d'acquisition et d'appariement d'un couple stereoscopique d'images comprenant 
un systeme d'acquisition simplifle permettant I'acquisition d'un couple pour un 
faible coefficient stereoscopique et un systeme de traltement mettant en oeuvre 
la methode d'appariement d'un couple stereoscopique d'images selon la 
presente invention. 
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PRESENTATION DE LMNVENTION 

invention concerne un ensemble d'acquisition et d'appariement d'un 
couple stereoscopique d'images comprenant un systeme d'acquisition d'un 
couple stereoscopique d'images et un systeme de traitement du couple 
stereoscopique acquis, caracterise en ce que le systeme d'acquisition du couple 
stereoscopique comporte un instrument d'acquisition unique comprenant deux 
capteurs CCD dans ie plan focal optique, chaque capteur CCD permettant 
Tacquisition d'une image du couple stereoscopique damages. 

Des aspects preferes, mais non limitatifs de i'ensemble d'acquisition et 
d'appariement d'un couple stereoscopique d'images selon I'invention sont les 
suivants : 

chaque capteur CCD est une barrette CCD ; 
chaque capteur CCD est une matrice CCD ; 
Ie systeme de traitement comporte : 

- des moyens aptes a traiter les deux images du couple stereoscopique, 
selon au moins une direction de traitement et selon differents niveaux de 
resolution du plus grossier au plus fin, lesdits moyens de traitement 
comprenant des moyens qui, pour chaque changement de niveau de 
resolution, sont aptes a : 

- determiner, en chaque point de la premiere image, une fenetre 
de correlation optimale, 

- calculer, par correlation avec les fenetres de correlation 
determinees, les ecarts de position entre chaque point de la 
premiere image et son homologue radiometrique dans la 
deuxieme image. 
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- obtenir, pour le niveau de resolution courant, une carte des 
ecarts de position entre chaque point de la premiere image et son 
homologue radiometrique dans la deuxieme image, 

- effectuer un recalage barycentrique des points de ladite carte 
obtenue. 

- des moyens aptes a obtenir, pour le niveau de resolution le plus fln, !a 
carte des ecarts de position entre chaque point de la premiere image et 
son homologue radiometrique dans la deuxieme image ; 

les moyens aptes a traiter les deux images du couple comprennent en 
outre : 

- des moyens aptes a effectuer une convolution des images du couple 
par une fenetre de convolution, 

- des moyens aptes a effectuer un zoom de taille 2 des images du couple 
stereoscopique ; 

les moyens aptes a determiner la fenetre de correlation optimale en 
chaque point de la premiere image comprennent : 

- des moyens aptes a calculer une courbure de la premiere image pour 
chaque taille de fenetre de correlation, 

- des moyens aptes a calculer une courbure de I'image de bruit pour 
chaque taille de fenetre de correlation, 

- des moyens aptes a calculer, pour chaque taille de fenetre de 
correlation, le rapport (SNRc) des courbures precedentes, 

et des moyens aptes a choisir la plus petite taille de fenetre de correlation 
telle que ledit rapport soit superieur a un parametre d^exhaustivite saisi 
par Tutilisateur. 

les moyens aptes a traiter les deux images comprennent en outre des 
moyens aptes a effectuer une interpolation de la deuxieme image grace 
a une carte des ecarts de position obtenue pour un niveau de resolution 
precedent. 
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les moyens aptes a traiter les deux images comprennent en outre : 

- des moyens aptes a rejeter les points homologues issus de fenetres de 
correlation plus grandes et contenant des fenetres de correlation plus 
petites, 

- des moyens aptes a rejeter les points homologues de rayons de 
courbure verifiant une condition de rejet apres interpolation de la 
deuxieme image, 

- des moyens aptes a completer iterative m en t, dans la carte des ecarts 
de position correspondant au niveau de resolution courant, des points 
d'ecarts de position lacunaires en fonction des points d'ecarts de position 
existants. 

- des moyens aptes a lisser la carte des ecarts de position correspondant 
au niveau de resolution courant par une fenetre de convolution. 

- des moyens aptes a composer la carle des ecarts de position calculee 
pour le niveau de resolution courant avec la carte des ecarts de position 
calculee pour le precedent niveau de resolution afin de mettre a jour 
ladite carte des ecarts de position calculee pour le niveau de resolution 
precedent. 

les moyens aptes a rejeter les points homologues verifiant une condition 
de rejet comprennent : 

- des moyens qui, pour le point courant, sont aptes a calculer sur la 
premiere image les courbures de ses deux points voisins suivant chaque 
direction de traitement, 

- des moyens qui, pour le point courant, sont aptes a calculer Tecart 
inferieur des courbures des 4 connexes au point courant de la premiere 
image, 

- des moyens qui, pour le point courant, sont aptes a calculer, dans la 
deuxieme image» son point homologue interpole par la carte des ecarts 
de position, 
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- des moyens qui, pour le point courant, sent aptes a oalculer I'ecart entre 
sa courbure et celle du point homologue, 

- des nnoyens qui, pour le point courant, sent aptes a rejeter le point 
homologue si cet ecart est plus important que I'ecart inferieur calcule sur 
la premiere image. 

L'invention concerne egalement un procede d'acquisition et 
d'appariement d'un couple stereoscopique dimages, caracterise en ce qu'il 
comprend [^acquisition d'un couple stereoscopique damages par un instrument 
d'acquisition unique comprenant deux capteurs CCD situes dans !e plan focal 
optique, le premier capteur permettant Tacquisition d'une premiere image "du 
couple, et le deuxieme capteur permettant ['acquisition d'une deuxieme image 
du couple. 

Des aspect preferes mais non limitatifs du procede selon la presente 
invention sont les suivant : 

le procede comporte en outre les phases consistant a : 

- traiter grace a des moyens de traitement les deux Images du couple 
stereoscopique, selon au moins une direction de traitement et selon 
differents niveaux de resolution du plus grossier au plus fin, en, pour 
chaque changement de niveau de resolution : 

- determinant, en cliaque point de la premiere image, une fenetre 
de correlation optimale, 

- calculant, par correlation avec les fenetres de correlation 
determinees, les ecarts de position entre chaque point de la 
premiere image et son homologue radiometrique dans la 
deuxieme image, 

- obtenant, pour le niveau de resolution courant, une carte des 
ecarts de position entre chaque point de la premiere image et son 
homologue radiometrique dans la deuxieme image, 
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- effectuant un recalage barycentrique des points de ladite carte 
obtenue. 

- obtenir, pour le niveau de resolution le plus fin, la carte des ecarts de 
position entre chaque point de la premiere image et son homologue 
radiometrique dans [a deuxieme image ; 

le precede comporte une phase consistant a traiter grace a des moyens 
de traitement les deux images du couple stereoscopique selon une 
direction epipolaire ; 

le procede comporte une phase consistant a traiter grace a des moyens 
de traitement les deux images du couple stereoscopique selon une 
direction epipolaire et une direction orthogonale a la direction epipolaire ; 
le procede comprend en outre une phase de saisie par un utilisateur d'au 
moins un parametre d*exhaustivite correspondant au choix d'un rapport 
signal sur bruit de correlation suivant chaque direction de traitement ; 
le procede comporte en outre une phase d'acquisition de donnees 
comprenant Tacquisition : 

- d'une image de bruit par chaque capteur CCD, 

- du type de la fenetre de correlation, 

- d'une direction epipolaire pour chaque image du couple pour un 
coefficient stereoscopique moyen, 

- d'au moins une valeur de Fecart de position maximum suivant chaque 
direction de traitement ; 

le procede comprend en outre une phase de pretraitement des donnees 
consistant a : 

- calculer un nombre de niveaux de resolution suivant chaque direction 
de traitement en fonction de. chaque valeur de Pecart de position 
maximum, 

- filtrer les images du couple stereoscopique d'images par convolution 
avec une fenetre de convolution, 
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- filtrer I'lmage de bruit par convolution avec ladite fenetre de convolution, 

- interpoler dans une geometrle epipolaire le couple stereoscopique 
d'images ; 

la phase de traitement des images comprend en outre les etapes 
consistant a : 

- effectuer une convolution des images du couple par une fenetre de 

convolution, 

- effectuer un zoom de taille 2 des Images du couple stereoscopique ; 
I'etape de determination de la fenetre de correlation optimale consiste a, 
pour chaque point de la premiere image : 

- calculer une courbure de la premiere image pour chaque taille de 
fenetre de correlation, 

- calculer une courbure de i'image de bruit pour chaque taille de fenetre 
de correlation, 

- calculer, pour chaque taille de fenetre de correlation, le rapport (SNRc) 
des courbures pr§cedentes, 

et a choisir la plus petite taille de fenetre de correlation telle que ledit 
rapport • soit superieur a chaque parametre d'exhaustivite saisi par 
I'utiiisateur. 

la phase de traitement des images comprend en outre une etape 
consistant a effectuer une interpolation de !a deuxieme image grace a 
une carte des ecarts de position obtenue pour un niveau de resolution 

precedent. 

la phase de traitement des images comprend en outre les etapes 
consistant a : 

- rejeter les points homologues issus de fenetres de correlation plus 
grandes et contenant des fenetres de correlation plus petlfes. 

- rejeter les points homologues de rayons de courbure verifiant une 
condition de rejet apres interpolation de la deuxieme image. 
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- completer ite ra tivement, dans la carte des ecarts de position 
correspondant au niveau de resolution courant, des points d'ecarts de 
position lacunaires en fonction des points d^ecarts de position existants. 

- lisser la carte des ecarts de position correspondant au niveau de 
resolution courant par une fenetre de convolution. 

- composer la carte des ecarts de position calculee pour le niveau de 
resolution courant avec la carte des ecarts de position calculee pour le 
precedent niveau de resolution afin de mettre a jour ladite carte des 
ecarts de position calculee pour le niveau de resolution precedent. 
I'etapede rejet des points homologues verifiant une condition de rejet 
consiste a, pour le point courant : 

- caiculer sur la premiere image les courbures de ses deux points voisins 
suivant chaque direction de traitement, 

- caiculer I'ecart inferleur des courbures des 4 connexes au point courant 
de la premiere image, 

- caiculer. dans la deuxieme image, son point homologue interpole par la 
carte des ecarts de position, 

- caiculer i'ecart entre sa courbure et celle du point homologue, 

- rejeter (e point homologue si cet ecart est plus important que Tecart 
inferleur calcule sur la premiere image. 

la phase d'obtention, pour le niveau de resolution le plus fin, de la carte 
des ecarts de position, comprend Tobtention : 

- d'au moins un tableau representant la carte des ecarts de position 
suivant chaque direction de traitement, 

- d'au moins un tableau representant une carte des tallies des fenetres de 
correlation employees en tout point non rejete de la premiere image, 
suivant chaque direction de traitement 

les fenetres de convolution et de correlation utilisees sont des fenetres de 
convolution et de correlation de type prolate. 
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PRESENTATION PES FIGURES 

D'autres caracteristiques et avantages de ['invention ressortiront encore 
de la description qui suit, laquelle est purement illustrative et non limitative et 
doit etre lue en regard des dessins annexes, sur lesquels : 

- La figure 1 illustre une scene representee dans un couple stereoscoplque 
d'images; 

- La figure 2 illustre une vue d'une sc§ne de la composante en colonnes, 
d'un decalage sinusoTdal injecte, et du resultat de la mesure de ce 
decalage par correlation ; 

- La figure 3 illustre une vue en perspective d'un systeme de stereoscopie ; 

- La figure 4 illustre un graphique representant le taux de correlation 
suivant une direction d'une image de type SPOTS en fonctlon de I'ecart 
de position entre le point courant de la premiere image et le point courant 
de la deuxieme Image, le point courant correspondant dans chaque 
image au point corr6l§ ou point au centre de la fenetre SPOTS (cf. p 23) ; 

- La figure 5 illustre un graphique representant le taux de correlation d'une 
fenetre de type hypomode carre 2X2 en fonction de I'ecart de position 
entre le point courant de la premiere image et le point courant de la 
deuxieme Image, le point courant correspondant dans chaque image au 
point correle ou point au centre de la fenetre ; 

- La figure 6 illustre un graphique representant le taux de correlation d'une 
fenetre de type prolate en fonction de I'ecart de position entre le point 
courant de la premiere image et le point courant de la deuxieme image, le 
point courant correspondant dans chaque image au point correle ou point 
au centre de la fenetre ; 

- La figure 7 illustre une vue d'une scene, des points de la scene 
conserves apres application d'un critere en ligne, des points de la scene 
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conserves aprds application d'un critere en colonne, et dee points de la 
scene conserves apres application d'un critere en iigne et en coionne ; 
La figure 8 est un grapliique de la variation de la precision aitimetrique en 
fonction du coefficient stereoscopique pour quatre metliodes differentes 
d'appariement d'images stereoscopiques (fenetres creneaux. prelates, 
prelates + correction barycentrique) ; 

La figure 9 illustre une vue d'un instrument dedie a la stereoscopie ; 

La figure 10 illustre une vue en perspective de deux images formant un 

couple stereoscopique d'images et de la scene observee ; 

La figure 1 1 illustre une vue d'un couple d'images presentant chacune 

trois reponses impulsionneiles ou fonction de transfert de modulation 

(FTM). 



DESCRIPTION D'UN MODE DE REALISATION DE L'lNVENTIQN 

La methode d'appariement fin d'images stereoscopiques selon la 
presente invention va maintenant §tre detalllee en reference aux figures 1 a 12. 
Les elements equivalents representes dans les differentes figures porteront les 
m§mes references numeriques. 

Le proced§ ci-presente permet le calcul de cartes precises des ecarts de 
position entre images de couples stereoscopiques ayant des coefficients 
stereoscopiques falbles avec la m§me precision aitimetrique que pour des 
coefficients stereoscopiques forts. Ce precede fonctionne jusqu'a des 
coefficients stereoscopiques tres faibles (0.01) sans degrader la precision 
aitimetrique. 

Cette configuration stereoscopique d'acquisition (b/h faible) presente 
I'avantage de limiter les parties cacfiees, c'est-a-dire des parties n'apparaissanf 
que dans I'une des deux images, et par exemple de voir le fond des rues en 
viile. En effet, comme I'illustre la figure 3. les parties cachees 14, sent liees a 
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Tacquisition par un capteur CCD 15 ou 16 d'une prise de vue oblique 2D d'un 
objet 3D. Aussi, i*acquisition sous coefficient stereoscopique faible permet une 
mise en oeuvre plus automatique de la realisation des cartes des ecarts de 
position entre couples damages en evitant la gestion, toujours laborieuse, des 
parties cachees. 

Las differentes etapes ayant conduit a la determination du precede et du 
systeme vont maintenant etre presentees. 

Les modeles utilises en traitement d'image sont des modeles discrets. 
Gependant, il est souvent plus facile de commencer par raisonner sur un 
modele continu, et de passer ensuite au cas discret. C'est cette approche qui 
est utilisee dans la suite pour decrire ia methode d'appariement d'image. 

1. Modelisation continue de la correlation unidirectionnelle 

La determination du modele continu de la correlation unidirectionnelle va 
conduire a Tobtention d'une equation permettant de relier la mesure de Tecart 
de position effectuee par I'operation de correlation suivant une direction de 
traitement en fonction de la valeur reelle de Tecart de position. Cette equation 
mettra en evidence le defaut precedemment indique de la correlation, dit 
"adherence". Un traitement permettant de llmiter ce phenomene d'adherence 
sera ensuite presente. Enfin, une seconde relation, a savoir une condition 
morphologique directionnelle locale dont depend la precision de Tecart de 
position mesure sera egalement obtenue a Tissue de ia modelisation continue 
de la correlation unidirectionnelle, 

Dans cette partie, il est fait Thypothese d'une direction unique des cartes 
des ecarts de position, par exemple suivant des droites dites epipolaires, et dont 
une definition va maintenant etre donnee. Comme illustre a la figure 10, dans un 
couple stereoscopique, on definit Pensemble des plans epipolaires comme etant 
I'ensemble des plans passant par les deux centres optiques CI et C2. 
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Les points P, P1 et P2 font respectivement partis de la scene 38, de la 
premiere image 36, et de la deuxieme image 37 du couple stereoscopique. II 
apparaTt que pour tous les points de la scene 38, le point P et ses points images 
P1 et P2 appartiennent au meme plan epipolaire 35. Les droites intersection des 
plans images 36 et 37 avec le plan epipolaire 35 sont les droites epipolaires. 
Celles-ci sont connues puisque les positions des capteurs CCD C1 et C2 sont 
connues. 

La formulation continue du coefficient de correlation non lineaire est 
effectuee dans le cas de couple stereoscopique d'Images faiblement bruitees. 
Le coefficient de correlation suivant la direction du vecteur unite v dans un 
voisinage centre en {t} est : 

p (^hi) - ' ^^^^ ' ^^^^^ ~ ^^^^ ^ ^ 

( - est I'operateur de correlation, 1 est Umage de reference, / est rimage a 
apparier et cp la fenetre ). 

On cherche Tecart de position vectoriel At qui maximise ce coefficient de 
correlation. En effet, Tecart de position pour lequel la correlation est maximale 
correspond a I'ecart de position entre les coordonnees d*un point de [a premiere 
image, et les coordonnees de son homologue radiometrique dans la deuxieme 
image, 

Soit u(t) la fonction vectorielle qui a tout point t de Timage associe At. Ce 
coefficient de correlation est calcule a partir d'un rapport, il est done non lineaire. 
Les ecarts de position etant des differences de positions locales relatives, la 
recherche de i'ecart de position maximisant le coefficient de correlation peut etre 
ramenee a une recherche autour du point to=0. 
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On obtlent : 

j(p-I ■ I(x + UQ)dx 



^ j^-I-(x)dxj^-I^(x + Uo)dx 

Mt. ®st une constante et uO = u(0) = At, I est I'image de reference, / est 
Timage a apparier et (p la fen§tre . 

L'objectif est de rechercher uO tel que ce coefficient de correlation soit 
maximunn. La formule I(x) = I(x+s(x)) permet le passage d'une image du 
couple stereoscopique a I'autre par un « petit » deplacement de pixels (on 
donnera un sens particulier au terme « petit » par la suite). Un pixel est 
Telement le plus petit d'une Image auquel on puisse associer indlviduellement 
un niveau de gris representatif d'une intensite. 

En effectuant un developpement limite au deuxieme ordre, suivant la 
direction du vecteur unitaire v centre en {0} et apres remplacement de / par sa 
valeur approchee en fonctlon de/, le coefficient de correlation devient : 

mi 

ou r est la derive de I suivant v . 
Equation reliant la mesure de I'ecart de position a sa valeur reelie 



La recherche du maximum du coefficient de correlation precedent conduit 
a rechercher I'annulation dans une direction quelconque de traitement de la 
derivee de ce coefficient. 
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Solf u la mesure de la valeur de i'ecart de position suivant la direction v 
au point t correspondant au maximum du coefficient de correlation (p'=0) et soit 
8 I'ecart de position reel suivant cette meme direction. 

L'equatlon fondamentale de la correlation non ilneaire est alors 

ou C^-(p*d est la courbure de correlation suivant la direction v, avec (p la 
fenetre de correlation et rf(/,7') = ^C_:^^ est la quasi-densite de correlation 

\cp ^ IP 1 

®*^(0=-— . Elle ne depend que de ('image de reference I et de sa derivee. 

II 11^^ 

Toutes les derivees sont exprimees dans la direction de traitement v . 

Cette relation fondamentale de la correlation permet de reiier la mesure 
de I'ecart de position effectue par I'operation de correlation a la valeur reelle de 
I'ecart de position, et ce sans passer par le caicul explicite de la correlation, 
calcul toujours tres couteux en temps. 

On remarque, en observant la formuie permettant d'obtenir la quasi- 
densite de correlation, que celle-ci est principalement . fonction de l'^. Cette 
quasi-densite de correlation est done tres sensible aux contours, et plus 
generalement aux singularites (points oCi certaines quantites deviennent 
infinies). Pour pallier le probleme d'adherence susclte, on utilise une methode 
dite de recalage barycentrique permettant de revenir a I'ecart de position vral. 
Cette methode va maintenant §tre decrite. 



Methode de recalage barycentrique 
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La methode dite de recaiage barycentrique prend sa source dans [a 
constatation sulvante. L'equation fondamentale au point de coordonnee, par 
example 0, devient si I'image a une quasi-densite concentree en {xo'}, 



Ceci signifie que la mesure au point 0 est egale a Tecart de position vrai 
au point {xo'} (au signe pres). II s'agit precisement du phenomene d'adherence 
cite plus haut. Le phenomene d'adherence de Umage affecte au point courant 
Tecart de position vrar fssu majoritairement d'un autre point plus singulier. La 
methode de recaiage barycentrique va done consister a rechercher la position la 
« plus singuliere », et a lui affecter la valeur de Tecart deposition mesuree. 

La fenetre de correlation q> est positive et Ton se place en 0. On 
recherche le barycentre sur le support de (p dont la quasi-densite est (p{x)d{x) . 
II s'agit de la position : 



ou le vecteur OPdecrit (e support de ^. Si la masse est concentree 
autour du centre de gravite aiors Tintegrale se restreint au voisinage du centre 
de gravite v(pG). De plus, I'ecart de position peut etre considere comme 
constant dans ce petit voisinage. 




OG = 




D'ou 




Done 
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Tous les termes de la relation precedents sont connus et done 
calculables. [| suffit d'affecter la valeur mesuree du decalage u(0) au point de 
coordonnees G. 

Une relation permettant de relier la mesure de I'ecart de position 
effectuee par I'operation de correlation a la valeur reelle de I'ecart de position 
ainsi qu'un traitement permettant de palier au probleme d'adlierence de la 
correlation viennent done d'etre presentes. Comme indlque dans la presentation 
du modele, une seconde relation, a savoir une condition morphologique 
directionnelle locale va maintenant etre determinee. 

Condition morphologique directionneUe locale 

La recherche du maximum de la surface de correlation permet de 
retrouver a partir d'un point de la premiere image son homologue radiometrique 
dans la deuxieme image. 

Cette recherche du maximum implique que la deriv^e seconde du 
coefficient de corr§lation soit strictement negative afin qu'il n'y ait pas plusieurs 
maximums a la surface de correlation. En effet, la presence de plusieurs 
maximums ^ la nappe de correlation rendrait ambigue la determination du point 
homologue. 

L'expression de cette derivee seconde suivant la direction v est appelee 
relation morphologique. 

/?"=C est appele courbure de correlation suivant la direction v . La courbure doit 
etre suffisante, et une expression quantitative en fonction du bruit va maintenant 
etre donnee au terme "suffisante". 
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En admettant que la nature du bruit soit connue, par example un bruit 
Gaussien d'ecart type a et que ce bruit ait ete regularise par une Gaussienne g 
(ecart type typique de 0,5 pas de quantification) 

^ -lo-Hgiljuig?) 

Q = rrr ' ^(^) rintegrale de 9. 



Alors le rapport signal a bruit de courbure de correlation se definSt par : 

C 



C 



qui est independant de I'image. 



La condition de correlation est : SNRc superieur a un seui! de i'ordre 
preferentlellement de 10 (un ordre de grandeur) qui permet de clioisir les points 
retenus. Si cette condition est rempHe, le bruit peut etre neglige dans le calcul 
des relations fondamentales et morphologiques, et le nnodele etabli plus haut est 
applicable. 

Une condition nnorphologique vient done d'etre etablle. Cette condition 
morphologique ne depend pas de Tecart de position. II s'agit d'un resultat d'une 
importance pratique non negligeable. En effet, il est possible de verifier en tout 
point de IMmage la realisation de cette condition et ceci sans recourir au calcul 
laborleux de la correlation elle-meme. Ce resultat justifie a posteriori le choix du 
coefficient de correlation pondere. 

On a done effectue ici la modelisation continue de la correlation non 
lineaire unidirectionnelle. 11 reste a passer au cas discret puisque une image est 
une suite de points et par consequent un modele discret. 



2. Modelisation discrete de la correlation non lineaire unidirectionnelle 
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La modelisation fine du calcul du coefficient de correlation discret conduit, 
comma il va etre precise ici, a : 

- effectuer un zoom de facteur 2 sur ie couple d'images, 

- choisir un type de fonction comme fenetre de correlation, 

- interpoler les images d'entree et la nappe de correlation avec un filtre 
sinus cardinal. 

Enfin un test d'arret des Iterations pour une direction de traitement et un 
calcul multi niveaux de resolution, dit multiechelle, seront effectues. 

Echantillonnaae 

Apres avoir effectue un echantillonnage de bonne qualite, c'est-a-dire 
*S'=:i?, ou S est Ie support spectral de Timage et R la cellule reclproque 
correspondant a Techantillonnage, on calcule Ie coefficient de correlation entre 
les images numerisees d'un couple stereoscopique. Ce calcul necessite de bien 
exprimer 11' et 1'^ qui "vivent" spectralement dans R+R que Ton nommera tout 

simplement 2R. Le calcul du coefficient de correlation Po(^o ) = ' ^^^^( ^J±L 



a partir de son numerateur et de son denominateur necessite un bon 
echantillonnage de I'image au carre. 

Une solution preferee pour calculer le coefficient de correlation est de 
calculer separement numerateur et denominateur, puis d'en realiser le quotient. 
En effet, le calcul direct fait intervenir les derivees au carre par exemple dans 
I'expression faisant intervenir IXx)^ (cf. relation fondamentale). 

Ces deux calculs ne necessitent bien evidemment qu'en apparence le 
meme support spectral 2S. Le calcul direct s'effectue sur des termes bien plus 
Haute Frequence (HF) (terme en 1'^) que celui ou denominateur et numerateur 
(terme seulement en P) sont calcules separement avant d'en faire le quotient. 
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De plus, le calcul du denominateur fait intervenir un carre de Timage. 

Pour satisfaire aux conditions d'ecliantillonnage precedemment decrites, 
un zoom de deux sur ies images du couple stereoscopique est realise. Ce zoom 
consiste, pour cliacune des images du couple stereoscopique, a doubier le 
nombre de points en ligne en inserant un point de valeur inconnue entre deux 
points de I'image de valeurs connues, a doubier le nombre de points en colonne 
en inserant un point de valeur Inconnue entre deux points de I'image de valeurs 
connues, et a determiner par interpolation Ies valeurs inconnues en utilisant Ies 
valeurs connues. En effet, I'interpolation sur le zoom de deux est une estimation 
de valeurs intermediaires dans une serie de valeurs connues. L'interpolatlbn sur 
le zoom de deux des images doit etre precise. 

La m§thode utilise preferentlellement des filtres separables d'interpolatlon 
longs de type sinus cardinal pour Ies diverses interpolations. La taille du filtre 
sinus cardinal d'interpolation de I'image pour effectuer un zoom sera 
preferentlellement de tailJe 35. Un filtre de type sinus cardinal sera egalement 
utilise pour realiser Ies interpolations sur la nappe de correlation. On prendra 
preferentiellement un filtre sinus cardinal de taille 1 1 pour Ies interpolations sur 
la nappe de correlation lors de la recherche sous pixellaire du maximum de 
correlation. 

Une fols le zoom de deux realise, on calcule le coefficient de correlation 
entre Ies images numerisees du couple stereoscopique. Comme deorit 
precedemment, une fenetre 3 de correlation est necessaire pour calculer !e 
coefficient de correlation entre Ies deux images du couple. Les fenetres de 
correlation peuvent etre de differents types. Cependant, certains types de 
fenetres minimisent I'impact de I'erreur sur la mesure des ecarts de position. 
L'analyse ayant conduit a la determination d'un type prefere pour la fenetre de 
correlation va maintenant etre presentee. 

Determination d'un tvpe pr6fere pour la fenetre de correlation 
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Au voisinage du maximum, la nappe de correlation doit etre convexe afin 
d'assurer la convergence de Talgorithme de recherche du maximum principal. 
En effet, si la nappe de correlation n'est pas convexe, cela signifie que plusieurs 
points de la deuxieme image du couple d'images peuvent potentiellement 
correspondre au point courant dans la premiere image. 

Or un seui de ces points est le juste point, c'est-a-dire rhomologue 
radiometrique vrai, du point courant de la premiere image. Pour une image de 
reponse impulsionnelle connue, le voisinage de taille maximum assurant cette 
propriete de convexite du maximum de la nappe va done malntenant §tre 
recherche. 

L'element le plus fin qui puisse §tre trouv§ sur une image est egal a la 
reponse impulsionnelle. Des lors, il suffit d'etudier numeriquement la forme de la 
nappe de con-elation pour un couple d'images identiques reduites a trois 
reponses impulsionnelles distantes de la longueur Al appartenant a I'intervalle 
en pixels [Q;7], et de rechercher la distance maximum entre les deux reponses 
impulsionnelles pr6servant une nappe de correlation convexe. Mors, cette 
valeur correspond a ['exploration maximale possible et elle donne un sens a la 
notion d'ecart de position dit "petit" utilisee en introduction du modele d'images. 

L'analyse 1D est suffisante dans le cas de reponses impulsionnelles 
separates. Pour realiser l'analyse 1D, on prend deux images 1,2 identiques 
pr§sentant des reponses impulsionnelles 301, 302, 303 distantes de Al 304, 
comme illustre a la figure 12. Le fait de prendre deux images identiques pour 
former le couple stereoscopique d'images induit que les homologues 
geometrique et radiometrique dans la deuxieme image d'un point de la premiere 
image sont confondus. On fait ensuite glisser suivant la direction d'analyse, la 
fen§tre 3 de correlation le long de la deuxieme image 2. 
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Pour chaque deplacement de 1 point de la fenetre de correlation, on 
calcule la correlation. On obtient las graphiques des figures 4, 5 at 6 presentant 
le taux de correlation entre les deux images identiques en fonction de la position 
de la fenetre 3 de correlation suf la deuxieme image. 

L'abscisse x=0 correspond au cas 300 ou la fenetre de correlation est 
placee sur Thomologue geometrique du point courant de la premiere image, 
c'est-a-dire que I'on mesure le coefficient de correlation entre un point de la 
premiere image et son homologue radiometrique vrai dans la deuxieme image. 

L'abscisse x=1 correspond au cas ou la fenetre de correlation est 
decalee dans la deuxieme image de 1 point par rapport a la position de 
rhomologue radiometrique vrai. 

Uabscisse x=2 correspond au cas 310 oO la fenetre de correlation est 
decalee dans la deuxieme image de 2 points par rapport a la position de 
i'homologue radiometrique vrai, et ainsi de suite. 

Les resultats de la correlation sont donnes aux figures 4, 5 et 6 pour les 
cas SPOT 5 JHR (figure 4), Hypomode standard (agglomeration par paquet de 
2x2 des pixels tout en conservant le meme echantillonnage) (figure 5) et prolate 
(figure 6), 

Comme on peut le remarquer a la figure 4, avec une fenetre de 
correlation de type SPOTS, Tecart de position maximum entre deux reponses 
impulsionnelles pour conserver une nappe convexe est d'un pixel. II faudra done 
avoir une zone d'exploration ne depassant pas le pixel. En effet, pour la 
distance zl/ = 1, on a une nappe de correlation convexe. Par centre, pour la 
distance Al = 2, on a plusieurs maximums a la nappe de correlation. 

Comme illustre a la figure 12, ces maximums correspondent au cas 300 
ou la fenetre de correlation est placee sur Thomologue radiometrique vrai du 
point courant de la premiere image, et a d'autres cas oCi la fenetre de correlation 
est decalee par rapport a rhomologue radiometrique vrai. 
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On peut observer quatre maximums secondalres pour x=2, x=-2, x=4 et 
x=-4. Le cas 310 est obtenu pour x=2, c'est-a-dire lorsque ia fenetre de 
correlation est decalee par rapport a I'homologue vrai dans la deuxieme image 
2. Ce decalage est de deux points suivant la direction de traitement 

Le cas 320 est obtenu pour x=-2, c'est-a-dire lorsque dans la deuxieme 
image, ie decalage de la fenetre de correlation par rapport a I'homologue vrai 
est de deux points dans !a direction opposee a la direction de traitement. Les 
cas 320 et 330 sont obtenus pour x=4 et x=-4. 

Pour ces quatre cas 310. 320, 330 et 340, on a quatre maximums ^ la 
nappe de correlation. La nappe de correlation nous donne done cinq points 
homologues au point courant de la premiere image. Or un seul de ces 
homologues mesures est I'homologue vrai du point courant de la premiere 
image. C'est la raison pour laquelle 11 est necessaire d'assurer la convexite de la 
nappe de correlation. Avec une fenetre de type SPOT 5, I'ecart de position 
maximum pour conserver une nappe de correlation convexe est done de 1 pixel. 

Comme on peut le remarquer a la figure 5, I'hypomodlsation standard 
permet une exploration de 2 pixels. Cela signifie que I'ecart de position 
maximum que Ton peut observer entre deux points homologues est de deux 
pixels. 

Comme on peut le remarquer a la figure 6, la fonction prolate permet une 
exploration de 4 pixels. 

La solution prolate assure la mellleure convexite elargie de la nappe de 
correlation et de fortes proprietes de continuite sur la mesure des ecarts de 
position. La fonction prolate possede la propriete d'etre la fonction positive dont 
les supports sont les plus concentres simultanement en espace et en frequence. 
Une methode preferee utilisera done une fonction prolate comme fenetre de 
correlation. 

II vient d'etre demontre que certains types de fenetres minimisent I'impact 
de I'erreur sur la mesure des ecarts de position. Cette recherche de la precision 
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de mesure conduit a opter pour une fenetre de correlation egale a la fonction 
prolate. Par aifleurs, le phenomene d'adherence est d'autant plus faible que !es 
fenetres sent de support reduit. Pour minimiser le phenomene d'adherence, on 
utilisera done une fonction prolate de support le plus petit possible compatible 
avec les conditions d'application du modeie de correlation fine precedemment 
decrit et rappele ci-dessous. 

Determination en chaque point d^une taille optimale de la fenetre de correlation 

II a precedemment ete defini un seuil morphologique SNRc lie a la 
precision de la correlation, De plus, rimpossibllite de resoudre Tequation 
fondamentale a ete indiquee. Ceci conduit a rechercher des fenetres de support 
le plus petit possible compatible avec ce seuil SNRc. 

Cette approche conduit a rechercher la prolate de support spatial le plus 
faible, c'est-a-dire de support spectral le plus important, remplissant la condition 
morphologique dans la direction de traitement. Une adaptation de la fenetre en 
tout point pour remplir au mieux la relation morphologique fondamentale va 
maintenant etre proposee : 

- Soit la suite de fonctions prolate {Pn} associee a la suite de supports 
spectra ux circulaires decroissants {(pn}. 

- Soit cpQ egal a C le cercle inscrit dans 2R support spectral le plus 
grand possible sur Timage de zoom 2. 

Alors, la recherche de la fenetre optimale se resume a trouver en tout 
point la prolate de support spectral circulaire maximal appartenant a la suite 
{Pn} et verifiant la condition morphologique. Ceci est realise simplement par 
essais successifs des fonctions prolate de supports spectraux decroissants. 
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Traitement multiechelle 

La limite d'exploration imposee par le type de fenetre de correlation d'une 
part, et par la recherche en tout point de la taille de fenetre la plus petite d'autre 
part, ne permettent pas de determiner les points homologues radiometriques 
d'un couple stereoscopique lorsque Tecart de position entre ces points est 
important. 

Ce resultat conduit a envisager un traitement en niveaux de resolution 
(multiechelle), car il permet par variation de la resolution de {'image de preserver 
les conditions de convexite de la nappe quelle que soit Famplitude de Pecart de 
position recherche. Une variation de la resolution est une variation de la 
definition du pixel. 

Par exemple, si une image est representee sur un ecran de 21 cm en 
ligne et 18 cm en colonne. Et si cette image presents 210 pixels en iigne et 180 
pixels en colonne, alors cette image presente un niveau de resolution de 10 
pixels par centimetre, Une diminution du niveau de resolution (passage d'un 
niveau de resolution fin a un niveau de resolution plus grossler) va correspondre 
a une diminution du nombre de pixels en ligne et en colonne pour representer 
I'image. Pour une diminution de 10 du niveau de resolution, on n'aura plus que 
21 pixels en lignes et 18 en colonne, soit un niveau de resolution de un pixel par 
centimetre. 

L'approche multi niveau de resolution (multiechelle) dyadique est 
necessaire pour remplir la condition de convexite de la nappe de correlation 
quelle que soit Tamplitude des ecarts de position, Le couple d'images est 
degrade en resolution par convolution par une prolate de support spectral de 
rayon rc/s, ou « rc » est le rayon correspondant a Timage zoomee d'un facteur 2, 
et « s » est le niveau de resolution courant. 

La descents en niveau de resolution necessite de boucher les trous, ou 
points non renseignes des cartes des ecarts de position. Les trous sont bouches 
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iteratfvement par convorution avec une prolate circulaire cie rayon 7 par 
exemple. A cheque iteration les points manquants sont tour a tour renseignes. 
Ainsi est calcule une carte exhaustive des ecarts. L'iteration a chaque niveau de 
resolution consiste par interpolation successive a corriger geometriquement a 
['aide de cette carte suivant une direction, par exemple la direction epipolaire, !a 
deuxieme image pour la rendre de plus en plus semblable geometriquement a la 
premiere image. 

En outre, if est possible - si le cas epipolaire n'est pas pur - d'explorer 
successivement les directions orthogonales (ligne et colonne), de calculer 
successivement des ecarts de position monodlrectionnels par correlation, et 
d'effectuer les interpolations unidirectionnelles correspondantes. La 
connaissance des ecarts de position maximaux en ligne et colonne induisent 
une profondeur de traitement en niveau de resolution potentieilement differents 
suivant les llgnes et les colonnes. Le nombre de niveaux de resolution 
dyadiques necessaires est : 

- en ligne : si=log2(maxJ), 

- et en colonne : Sc=log2(max_c), 

o(i maxj et max_c sont les ecarts de position maximaux respectivement en 
ligne et colonne. 

A chaque niveau de resolution, les operations d'interpolation de la 
deuxieme image, suivant les cartes des ecarts de position monodirectionnels 
successives calcuiees, necessitent I'application d'une formule de composition 
des cartes des ecarts de position. Apres chaque iteration est effectuee la 
composition entre la carte des ecarts de position entre I'image interpolee et la 
premiere image et la carte des ecarts de position qui a permis de generer 
I'image interpolee (c'est-a-dire la carte des hearts de position de ['iteration 
pr^cedente). Ainsi, les erreurs d'interpolation ne se cumulent-elles pas a travers 
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les iterations car nous sommes toujours dans le cas de la composition de deux 
ecarts de position au maximum. 

La disparite composee avec recalage barycentrique apres deux iterations 
s'exprime de la maniere suivante (saciiant que la premiere est toujours 
regularisee et sans trous de correlation) : 

- Soit si la premiere carte des ecarts de position entre I'lmage 7 et la 
premiere image i. Cette fonction est continue et regulidre. 

- Soit 7i(n)=7 (n + sj(7j-^i)) I'lmage 7 corrlgee des ecarts de position 
de barycentre ou wj («) ^ -s^ (« -^^)ei 1 2 (m) = (m + S2im- ^2 )) , I'image 
Ii corrlgee des ecarts de position S2 de barycentre ^2 oCi 
M2 (/«) « -S2 (m-42)' Ml (») et U2 (m) sont les mesures d'ecarts de position 
mesures barycentriquement, alors : 

- hO^) = h(f'^ + ^20n-^2)=7 {m + S2(jm-42) + Sl(rn-\-£20n-'S2)-4l)) ■ 

Critere de relet des points mal apparies 

11 a ete precedemment fait mention de points non renselgnes. Un point 
non renseigne est un point de la premiere image dont I'appariement avec un 
point de la deuxidme Image a ete rejete. Le critere de rejet des points de fausse 
correlation consiste a comparer point a point les courbures de correlation entre 
ia premiere image et la deuxieme Image corrlgee des ecarts de position 
(interpolee dans la grille de reference). 

Un point est rejete a chaque niveau de resolution quand I'ecart entre les 
valeurs de courbures entres les deux images est superieur a une valeur. Cette 
valeur est egale au minimum d'ecart aux 4 points connexes du point courant de 
la premiere Image. 
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Enfin. pour chaque direction de traitement et pour le niveau de resolution 
de rimage, les points calcules a partir de prolate de supports contenant une 
prolate de supports plus petits sont elimines car la correlation en ces points est 
de precision obligatoirement plus faible. 

3. iVIethode de correlation fine bidirectionnelle 

Le cas de la methode de correlation fine unidirectionnelle a ete expose ci- 
dessus pour faciliter la comprehension du modele. Le cas de la methode de 
correlation fine bidirectionnelle va maintenant etre presents. 

Les images du couple stereoscopique sont re6chantillonn§es en 
geometrie eptpolaire suivant les lignes ou fes colonnes par interpolation. La 
recherche des tableaux de disparite 2D est effectuee alternativement suivant les 
lignes et les colonnes. La difficult^ majeure de restitution de la table des ecarts 
de position est generalement celle suivant la direction epipolalre. 

Les ecarts de position provoques par le relief sont en general tres haute 
frequence alors que ceux suivant la direction orthogonale sont souvent 
provoques par des vibrations du porteur et sont basse frequence. Ceci autorise 
des fen§tres de traitement plus grandes suivant la direction perpendiculaire aux 
epipolaires. Cette difference sera naturellement geree par le choix d'un rapport 
signal a bruit plus important suivant la direction orthogonale aux epipolaires que 
suivant la direction epipolaire. 

La condition fondamentale unidirectionnelle si v est la direction 
epipolaire est la relation unique : 

C^/Cj,>SNR^ 

La figure 7 pr§sente le masque 23 des points conserves apres 
application du critere de la courbure geometrique sur I'image de type campagne 
20 avec une fenetre 9x9 de type prolate. Ce masque 23 des points conserves 
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est obtenu en faisant rintersection du masque 21 des points conserves apres 
application du critere en ligne, avec le masque 22 des points conserves apres 
application du critere en colonne. 

Les zones noires de Timage correspondent aux points ayant ete rejetes 
apres application de la condition morphologique. Les points en blanc 
representent les points conserves. Les images de la figure 7 sent representees 
en noir et blanc pour une meilleure visualisation des points conserves. 

Comme I'illustre la figure 7, le masque 21 applique en ligne presente 
57.29% de points conserves, le masque 22 applique en colonne presente 
53,52% de points conserves, enfin, le masque 23 Intersection des deux 
precedents presente seulement 45.64% de points conserves, Cette Illustration 
confirme Tavantage non negligeable de Tutilisation du repere epipoJaire qui 
implique une condition morphologique unidirectionneile. 

4, Methode de correlation fine multidirectionnelle 

Uimage ainsi interpolee est ensuite analysee suivant sa direction 
orthogonale. Ce procede est itere. 

5. Precede d'appariement mis en place a partir du modele de correlation fine 

Le procede de calcul base sur la methode de correlation fine 
precedemment decrite va maintenant etre detaille en reference a la figure 11 . Le 
procede de correlation s'adapte automatiquement a I'image des que Ton a fixe 
le rapport signal a bruit de courbure de correlation SNRc. 

Ce procede de calcul est mulfi niveaux de resolution (multiechelle) 
dyadique. Cela signifie que le traitement aura lieu pour differents niveaux de 
resolufion, Le nombre de niveaux de resolution est connu des que Tecart de 
position maxlmale le long de la ligne et le long de la colonne Test. 
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Le procede de calcul base sur la methode de correlation fine est 
multidirectionnel. Las tallies de fenetres de correlation (prolate) sont calcuiees 
automatiquement en tout point et a tous les niveaux de resolution avec des 
fenetres de correlation qui peuvent etre tres differentes dans la direction 
epipolaire et la direction orthogonale, Elles sont fixees par le rapport signal a 
bruit de courbure de correlation different suivant la direction epipofaire et la 
direction orthogonale. On procedera en niveaux de resolution croissants, c'est- 
a-dire du plus grossier au plus fin. 

A Donnees : 

La premiere etape du procede consiste a racquisition des donnees d'entree. 

Pour utiliser le procede, les donnees necessaires en entree sont les suivantes : 
" le couple stereoscopique d'images, a savoir la premiere image et la 
deuxieme image, Ces deux images sont acquises a Taide d*un systeme 
d'acquisition particulier decrit plus precisement dans la suite du 
document Ce systeme permet Tacquisition d'un couple damages pour un 
faible coefficient stereoscopique. Ce systeme d'acquisition est un 
instrument unique comprenant 2 capteurs CCD, chaque capteur 
permettant Tacqulsition d'une image du couple stereoscopique. 

- recart type du bruit et le type de fonction FTM (donnee) : a savoir 
« prolate », « spots », « hypomode »... 

- les images du couples sont exprimees en direction epipolaire ligne ou 
colonne, 

- la valour de Tecart de position maximum suivant les epipolaires et la 
valeur de I'ecart de position maximum suivant la diYection orthogonale. 
Ces valeurs representent les seuils maximums du decalage, suivant les 
deux directions de traitement, entre un point de la premiere image et son 
homoiogue dans la deuxieme image. Ces valeurs vont permettre de 
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determiner le nombre de niveaux de resolution necessaires pour trouver 
rhomologue dans la deuxieme image d'un point de la premiere image. 

B Parametres utiiisateurs 

A l^etape suivante, rutilisateur est invite a entrer un parametre 
d'exhaustivite pour le procede: 

- le choix d'un rapport signal sur bruit (SNR) talon de correlation suivant la 
direction epipolaire et le choix d'un SNR talon de correlation suivant la 
direction orthogonale : SNR_e, SNR_o. Preferentiellement on proposera 
SNR_e=10 et SNR_o=100. Ces parametres SNR_e et SNR_o vont servir 
de valeurs seuils dans la determination en chaque point de la premiere 
image, d*une taille optimale pour la fenetre de correlation suivant chaque 
direction de traitement. 

C Parametres alqorithmiques 

Un certain nombre d'autres parametres sont pre regies : 

- le rayon maximum de la prolate circulaire de correlation. Ce rayon est le 
seuil definissant la taille maximum autorisee pour la fenetre de 
correlation. Ce rayon maximum sera preferentiellement egal a 25 pour un 
zoom image de 2. 

- le nombre d'iterations pour chaque direction de traitement, 
preferentiellement egal a 2. 

- le nombre d'iterations, lors du comblement de parties lacunaires de la 
carte des ecarts de position. En effet, le procede permet de definir les 
ecarts de position entre les points de la premiere image et leurs 
homologues respectifs dans la deuxieme image, Tous les ecarts de 
position mesures sont places sur une carte des ecarts de position. 
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Cependant, parmi ces ecarts, ceux qui ne repondent pas a un critere de 
ressemblance sont rejetes. La carte des ecarts de position presents done 
des points pour (esquels I'ecart de position entre un point de la premiere 
image et son homologue dans la deuxieme image n'est pas renseigne. 
Ces points non renseignes de la carte des ecarts de position sont 
appeles parties facunaires, et le nombre d'iterations lors du comblement 
de ces points est preferentiellement egal a 10. 

- le zoom de 2 des Images du couple a tous les niveaux de resolution. 

D Initialisation 

A partir des donnees d'entree, certains parametres sont initialises. 
On effectue done le calcul du nombre de niveaux de resolution suivant la 
direction epipolaire S_e et la direction orthogonale S_o en fonction des valeurs 
des ecarts de position maximums (suivant les deux directions de traitement) 
obtenues lors de la phase A. Le nombre de niveaux de resolution dyadiques 
necessaire est calcule a Taide de la formule du logarithme a base 2 des ecarts 
de position maximums en ligne et colonne precedemment citee : 

- s^e=fog2(valeur maximum du decalage selon la direction epipolaire), 

- s_^o=log2(valeur maximum du decalage selon la direction orthogonale), 

Par ailfeurs, les images d'entree doivent etre filtrees, lors du premier 
passage dans la boucle principale de traitement du procede de correlation. Ceci 
permet d'attenuer le bruit des images d'entree. On effectue done un : 

- fiitrage du couple stereoscopique damages par convolution avec une 
fenetre de convolution. Gette fenetre de convolution sera 
preferentiellement une prolate circulaire de rayon egal a 7. Gette 
operation de fiitrage de la premiere et de la deuxieme image du couple 
permet de regulariser les points de ces images. 
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- filtrage cie Timage de bruit par convolution avec une fenetre de 
convolution. Cette fenetre de convolution sera preferentiellement de type 
prolate circulaire de rayon preferentiellement egal a 7. L'interet de cette 
operation est la encore de regulariser les points de Timage de bruit. 

- interpolation dans la geometrie epipolaire du couple d'images. Les 
images du couple stereoscopique sont reechantilionnees en geometrie 
epipolaire sulvant les lignes ou les colonnes par interpolation. 

E Procede pour chaque niveau de resolution 

On effectue alors le traitement des Images du couple stereoscopique en 
niveaux de resolution dyadique, alternativement si besoin, suivant la direction 
epipolaire et orthogonale a la direction epipolaire (fonction de la profondeur de 
traitement en niveaux de resolution suivant la ligne et la colonne). Pour le 
traitement, on commence par le niveau de resolution le plus* grossier et on 
passe progressivement au niveau de resolution le plus fin. Les etapes qui 
suivent seront effectuees pour chaque niveau de resolution et se divisent en 
deux categories qui sont la preparation des donnees et la correlation des 
images stereoscopiques : 

E.1 Preparation des donnees 

L'etape de preparation des donnees est effectuee en realisant les etapes 
suivantes : 

- la convolution des images du couple par la prolate (prolate circulaire de 
taille preferentiellement egale a 7 pour le niveau de resolution 1, prolate 
separable pour les autres). Cette convolution permet la degradation en 
resolution des images de depart du couple. 
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le zoom de deux de la premiere image, de la deuxieme image et de 
I'image de bruit. Ce zoom de deux est realise par interpolation a I'aide 
d'une fonotion de type sinus cardinale. L'operatlon d 'interpolation 
consiste a, pour chaque image interpolee :" 

o doubler le nombre de points en ligne en inserant un point de valeur 

inconnue entre deux points de I'image de valeurs connues. 
o doubler le nombre de points en colonne en inserant un point de 

valeur inconnue entre deux points de I'image de valeurs connues. 
o passer une fenetre ponderee de type sinus cardinale en chaque 
point de valeur inconnue et estimer la valeur de ce point gr§ce aux 
points de valeurs connues. 
la determination de la fen§tre de correlation optimale. Pour cela, on 
effectue le calcul en tout point de la premiere image du rapport signal sur 
bruit de courbure de correlation SNRc en fonction de la fenetre de 
correlation (prolate circulaire). Puis on choisit la plus petite fenetre de 
correlation, suivant la direction epipolaire et perpendiculaire, verifiant la 
condition C^/Ci,>sm^. Cette determination est effectuee en realisant, 
en chaque point de I'image, les etapes suivantes; 

o calcul de la courbure de correlation du bruit Cb sur I'image de bruit 
pour chaque taille de fenetre de correlation de type prolate. Ce 
calcul est realise en utilisant la formule de la courbure de 
correlation citee lors de la presentation du modele de correlation 
fine unidirectionnelle, 
o calcul de la courbure de correlation de la premiere image pour 

chaque taille de fenetre de correlation de type prolate, 
o calcul pour chaque fenetre de correlation de type prolate du 
rapport signal sur bruit de courbure de correlation (SNRc) rapport 
des courbures precedentes, 
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o chofx de la plus petite taille de fenetre de correlation de type 
prolate tel que \CJC^\ > SNR_^ et \CJC^\ > SNR_^ . 

E.2 Correlation des images stereoscopiques 

Une fois les donnees preparees, on entre dans I'etape de traitement des 
donnees. Cette etape comprend les sous etapes qui suivent Ces sous etapes 
sont reiterees pour chaque point de (Image et suivant chaque direction de 
traitement : 

Pour chaque point de la premiere image, et pour chaque direction de traitement, 
on effectue : 

- rinterpolation de la deuxieme image avec la carte de disparite 
cumulee trouvee a Piteration precedente, c'est-a-dire pour le niveau 
de resolution precedent. Cette interpolation de Timage 2 est faite a 
I'aide de la carte de disparite. 

ie calcul de Tecart de position par correlation dans la fenetre de 
correlation, entre Ie point courant de la premiere image et son 
homologue suppose dans la deuxieme image, Ce calcul est realise 
pour tous les points de ia premiere image et permet I'obtention d'une 
carte des ecarts de position intermediaire (c'est-a-dire une carte ou la 
technique de correction barycentrique notamment n'a pas ete 
reaiisee). 

- Ie recalage barycentrique a chaque niveau de resolution au point pres 
de la carte des ecarts de position (intermediaire). Ceci permet de 
limiter Timpact du phenomene d'adherence de la correlation 
precedemment presente. En un point de la carte des ecart de 
position, I'etape de recalage barycentrique consiste a determiner Ie 
barycentre des points contenus dans la fenetre de correlation ayant 
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servi pour la determination de I'ecart de position antra le point courant 
de la premiere Image at son homologue dans la deuxieme image. 
Puis a assigner a ce barycentre la valeur de I'ecart de position du 
point courant dans la carte des ecarts de position. 

- Pour une meme resolution : 

- le rejet des points issus de fenetres plus importantes et contenant des 
fenetres plus petites. En effet, la correlation en ce point a I'aide de la 
fenetre de correlation de plus grande taille est de precision 
obligatoirement plus faible. 

- le rejet des points homologues de rayons de courbure trop differents 
apres Tinterpolation de la position de la deuxieme image. Pour cela, 
on va effectuer pour chaque resolution et direction : 

o le calcul de ia courbure de correlation du point courant (point 

etudie) de la premiere Image, 
o le calcul de la courbure de correlation de I'homologue dans la 

deuxieme Image du point courant de la premiere image, 
o le calcul de la difference entre les deux courbures 

precedemment calculees (ecart entre les courbures de 

correlation du point courant et de son homologue dans la 

deuxieme image), 
o le calcul de I'ecart de courbure minimum entre les courbures de 

correlation des quatre voisins (dans la premiere image) au 

point courant et ia courbure de correlation du point courant. 
o le rejet du point liomologue si I'ecart de courbure entre le point 

courant et son homologue est plus important que I'ecart 

inferieur calcule sur les quatre volslns du point courant de la 

premiere image. 

le rebouchage iteratif (preferentiellement 10 iterations) des cartes des 
ecarts de position iacunaires en passant par les points existants. En effet. 
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on a vu que les points homologues d'ecarts de courbure trop differents 
etaient rejetes. Par consequent certains points de la carte des ecarts de 
position vont etre non renseignes. On effectue done un rebouchage 
iteratif de la carte des ecarts de position. De plus la correction 
barycentrique produit des trous. 
Pour cela, on effectue : 

o la convolution par une fenetre de convolution de type prolate, 
o le rebouchage des points non renseignes de la carte des 
ecarts de position par convolution par une prolate. 

- le lissage de la carte des ecarts de position par une prolate separable 
correspondant au niveau de resolution courant (circulaire au niveau de 
resolution de Tinnage). Ce Hssage permet de debruiter la carte des ecarts 
de position. 

- a la connposition de la table des ecarts de position courante avec celle 
calculee a la precedente iteration (cumul). Cette composition permet de 
mettre a jour la table des ecarts de position calculee pour le niveau de 
resolution precedent. 

On effectue le traitement pour tous les niveaux de resolution en eliminant 
les points calcules avec une fenetre de correlation prolate cdntenant une fenetre 
de correlation prolate plus petite et ceci suivant les deux directions de 
traitement. On obtient alors en sortie du precede la carte des ecarts de position 
. pour le niveau de resolution le plus fin. 

F Sorties 

On obtient a la derniere etapes du precede les donnees de sortie du precede, II 
s'agit de tableaux de tallies egales a la taille de la premiere et de la deuxieme 
image. Ges tableaux sont : 

- la carte des ecarts de position (2 tableaux), 
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- la carte des tallies de prelates (2 tableaux) en tout point renseigne (non 

rejete) de la premiere Image au niveau de resolution le plus grossier. 

Le precede presente ci-dessus est mis en oeuvre dans un ensemble 
d'acquisition et d'appariement d'un couple stereoscopique d'images. Get 
ensemble d'acquisition et d'appariement comprend un systeme d'acquisition et 
un systeme de traitement. 

Le systeme d'acquisition permet I'acquisition des donnees d'entrees du 
precede (etape A du procede). Ce systeme d'acquisition est par exemple un 
satellite et comprend des moyens de communication (pour la communication 
avec le systeme de traitement), des moyens de traitement (de type processeur), 
des moyens memoire (pour stocI<er les images acqulses) et un capteur 
optoelectronique (systeme optique + capteurs CCD). 

Le systeme de traitement permet quant S lul de realiser I'appariement du 
couple stereoscopique d'images. Le systeme de traitement est par exemple une 
station de travail comprenant des moyens memoires (RAM, ROM) connectes a 
des moyens de traitement tel qu'un processeur, des moyens de visualisation 
tels qu'un ecran d'affichage et des moyens de saisie tels qu'un clavier et une 
souris. Le systeme de traitement est connecte avec des moyens de 
communication afin de recevoir les images a apparier acqulses par le systeme 
d'acquisition. 

Une comparaison de la methode de correlation fine avec des metliodes 
ne faisant pas intervenir le zoom de 2 du couple stereoscopique d'images, la 
modelisation de la correlation avec la recalage barycentrique, et le calcul multi 
niveaux de resolution (ou multiechelle) du maximum de la nappe de correlation 
est proposee, a partir de simulations representees a la figure 7. 

Les resultats obtenus pour quatre methodes d'appariement d'images 
stereoscopiques sent presentes (fig. 8 precision en fonction du b/h). ces 
methodes etant : la methode de correlation standard par une fenetre de 
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correlation de type prolate 26, la methode de correlation standard par une 
fenetre de correlation de type constante 25, la methode de correlation fine par 
une fenetre de correlation de type prolate 27, et la methode de correlation fine 
par une fenetre de correlation de type prolate en enlevant les points issus de 
fenetres incluant des fenetres plus petites 28. 

Les images du couple sont des images de Marseille de pas 
d'echantillonnage un metre avec une fonction de transfert respectant quasfment 
le principe de Shannon (valeur proche de 0 a la frequence de coupure). Le 
rapport signal a bruit des images est egal a 100. 

L'acquisition des images du couple stereoscopique est matricielle, Ces 
images sont calculees a partir d'une orthophotographie echantillonnee a un 
metre et d'un modele numerique de terrain de precision submetrique couvrant la 
meme zone. 

Des couples d'images stereoscopiques sont generes pour plusreurs b/h. 
Ce coefficient prend des valeurs entre 0.005 et 0.25. Ce coefficient est provoque 
par le seul depointage en tangage. 

Comme represents a la figure 8, la precision avec la methode de 
correlation fine est deux fois meiileure que cede mesuree avec les methodes 
classiques. De plus, on peut considerer que la correlation fine est de precision 
constante pour des b/h entre 0.01 et 0.15. La precision altimetrique dont Tecart 
type est ici presente en ordonnee avec comme unite ie pixel est meiileure que le 
pixel. Cette methode conserve une precision subpixellaire avec de falbles 
coefficients stereoscopiques (b/h). 

Par ailleurs, la methode preconisee dite de correlation fine sans inclusion 
de fenetres permet de rejeter les points d'ombres qui sont souvent aberrants. Le 
taux de correlation est proche de 1 si Ton neglige le phenomene d'ombre. 

Le precede d'appariement damages stereoscopiques decrit ci-dessus 
permet done le traitement d'images stereoscopiques pour un faible coefficient 
stereoscopique (b/h faible) avec la meme precision altimetrique que pour des 
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coefficients stereoscopiques forts. Ce precede fonctionne jusqu'a des 
coefficients stereoscopiques tres faibles (0.01) sans degrader la precision 
altimetrique, 

Ce precede de correlation fine est mis en osuvre dans le systeme de 
traitement de I'ensemble d'acquisition et d'appariement d'un couple 
stereoscopique d'images. 

Ce precede permet une nouvelie conception des systemes spatiaux 
d'acquisition et de photographie aptes a acquerir le couple stereoscopique 
d'images. Un coefficient stereoscopique (b/h) proche de 0.02 limite les 
directions de visees tiomologues a une valeur inferieure a ±un degr§. 

Un nouveau type de systdme spatial d'acquisition permettant I'acqulsitlon 
du couple stereoscopique d'images est repr§sente a la figure 9. 

Contrairement aux systdmes d'acquisition de I'art anterieur, le systeme 
d'acquisition selon la presente invention comporte un instrument d'acquisition 
unique (non represente) comprenant un systeme optique unique (non 
represente) et deux capteurs CCD 31, 32 symetriques dans le plan focal 
optique. 

Chaque capteur CCD represente a la figure 9 est une barrette constituee 
de detecteurs 33a. 33b, 33c. Ces detecteurs 33a, 33b, 33c sent par exempie 
des phofodiodes electroniques sensibles CCD a la lumiere et capable de 
convertir le signal lumineux en un courant electrique proportionnel a I'intensite 
de cette lumiere. Ces detecteurs 33a, 33b, 33c sont places cote a cote, en ligne 
et ferment la barrette 33. Chaque detecteur 33a, 33b, 33c est charge de 
I'observ^ation d'un pixel. Chaque detecteur 33a, 33b. 33c capte la lumiere 
provenant d'un pixel de terrain. 

Le Principe de fonctionnement d'un tei systeme d'acquisition d'image 
stereoscopique dans le cas spatial est le sulvant. 
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Chaque barrette 33, 34 permet Tacquisition d'une ligne dimage. La 
barrette 33 permet Tacquisition d'une ligne de la premiere image, et la barrette 
34 I'acquisltion d'une ligne de la deuxieme image. 

La premiere image et la deuxieme image du couple stereoscopique sont 
done acquises ligne par ligne par les barrettes CCD 33, 34 au fur et a mesure 
de Tavancement du satellite sur son orbite (entre le temps t et t+At). 

Les lignes d'images acquises par les premieres et deuxiemes barrettes 
33, 34 sont stockees dans des moyens memoire. 

Une fois les deux images 1, 2 acquises, le systeme d'acquisltion envoie 
(par des moyens de communication sans fil) le couple stereoscopique damages 
1 , 2 au systeme de traitement, Ce systeme de traitement est preferentiellement 
base au sol et permet Tappariement des points du couple. 

Dans le cas d'un systeme d'acquisition spatiale les deux barrettes 
Introduites dans le plan focal peuvent etre remplacees par deux matrices 
(chaque matrice etant constituee d*un ensemble de detecteurs CCD repartis en 
lignes et colonnes) egalement situees dans le plan focal. 

Ceci peut etre realise egalement avec un appareil photographique non 
spatial. Dans le cas d'un appareil photographique, le systeme d'acquisition du 
couple stereoscopique d'images comportera un capteur optoelectronique 
comprenant un systeme optique unique et deux capteurs CCD symetriques 
dans le plan focal. 

Le precede d'appariement fin d'un couple stereoscopique d'images mis 
en cBuvre dans le systeme de traitement reste vaiable pour les systemes 
d'acquisition comprenant un instrument spatial unique constitue de deux 
matrices ou de deux barrettes des lors que les perturbations d'attitudes sont 
corrigees ou negligeables. 
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REVENDICATIONS 

1. Ensemble d'acqulsition at d'appariement d'un couple stereoscopique 
d'images (1, 2) comprenant un systeme d'acquisition d'un couple 
stereoscopique d'images et un systeme de traitement du couple 
stereoscopique acquis, 
caracterise en ce que 

le systeme d'acquisition du couple stereoscopique comporte un Instrument 
d'acquisition unique (30) comprenant deux capteurs CCD (31, 32) dans le 
plan focal optique, chaque capteur CCD (31, 32) permettant I'acquisltion 
d'une image (1 , 2) du couple stereoscopique d'images. 

2- Ensemble selon la revendlcation 1. caract6rise en ce que chaque capteur 
CCD est une barrette CCD (33, 34). 

3. Ensemble selon la revendlcation 1, caracterise en ce que chaque capteur 
CCD est une matrice CCD (33, 34). 

4. Ensemble selon la revendlcation 1 ou 2 ou 3, caracterise en ce que le 
systeme de traitement comporte : 

- des moyens aptes a traiter les deux images (1, 2) du couple 
stereoscopique, selon au moins une direction de traitement et selon 
differents niveaux de resolution du plus grossier au plus fin, lesdits 
moyens de traitement comprenant des moyens qui, pour chaque 
changement de niveau de resolution, sont aptes a : 

- determiner, en chaque point de la premiere image (1), une 
fenetre de correlation (3) optimale, 

- calculer, par correlation avec les fenetres de correlation 
determinees, les ecarts de position entre chaque point (40) de la 
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premiere image (1) et son homologue radiometrique (41) dans la 
deuxieme image (2), 

- obtenir, pour le niveau de resolution courant, une carte des 
ecarts de position entre chaque point (40) de la premiere image et 
son homologue radiometrique (41) dans la deuxieme image, 

- effectuer un recalage barycentrique des points de ladite carte 
obtenue, 

- des moyens aptes a obtenir, pour le niveau de resolution le plus fin, la 
carte des ecarts de position entre chaque point de la premiere image (1) 
et son homologue radiometrique dans la deuxieme Image (2). 

Ensemble selon Tune des revendlcations 4. caracterise en ce que les 
moyens aptes a traiter les deux Images du couple comprennent en outre : 

- des moyens aptes a effectuer une convolution des images (1, 2) du 
couple par une fenetre de convolution, 

- des moyens aptes a effectuer un zoom de tallle 2 des images (1, 2) du 
couple stereoscopique. 

Ensemble selon la revendication 4 ou 5, caracterise en ce que les moyens 
aptes a determiner la fenetre de correlation optimale en chaque point de la 
premiere image comprennent : 

- des moyens aptes a calculer une courbure de la premiere image pour 
chaque taille de fenetre de correlation (3), 

- des moyens aptes a calculer une courbure de Timage de bruit pour 
chaque taille de fenetre de correlation (3), 

- des moyens aptes a calculer, pour chaque taille de fenetre de 
correlation, le rapport (SNRc) des courbures precedentes, 
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et des moyens aptes a choisir la plus petite taille de fenetre de correlation telle 
que ledit rapport soit superieur a un parametre d'exhaustivite saisi par 
I'utilisateur. 

7. Ensemble selon Tune quelconque des revendications 4 a 6, caracterise en 
ce que les moyens aptes a traiter les deux images comprennent en outre 
des moyens aptes a effectuer une interpolation de la deuxieme image (2) 
grace a une carte des ecarts de position obtenue pour un niveau de 
resolution precedent. 

8. Ensemble selon la revendication 7, caracterise en ce que les moyens aptes 
a traiter les deux images comprennent en outre : 

- des moyens aptes a rejeter les points homologues (40 , 41) issus de 
fenetres de correlation (3) plus grandes et contenant des fenetres de 
correlation (3) plus petites. 

- des moyens aptes a rejeter les points homologues (40, 41) de rayons 
de courbure verifiant une condition de rejet apres interpolation de la 
deuxieme image (2). 

- des moyens aptes a completer iterativement, dans !a carte des ecarts 
de position correspondant au niveau de resolution courant, des points 
d'ecarts de position lacunaires en fonction des points d'ecarts de position 
existants. 

- des moyens aptes a lisser la carte des ecarts de position correspondant 
au niveau de resolution courant par une fenetre de convolution. 

- des moyens aptes a composer la carte des ecarts de position calculee 
pour le niveau de resolution courant avec la carte des ecarts de position 
calculee pour le precedent niveau de resolution afin de mettre a jour 
ladite carte des ecarts de position calculee pour le niveau de resolution 
precedent. 
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9. Ensemble selon Tune des revendications 8, caracterise en ce que des 
moyens aptes a rejeter les points homologues verifiant une condition de rejet 
comprennent : 

- des moyens qui, pour le point courant, sont aptes a calculer sur la 
premiere image (1) les courbures de ses deux points voisins suivant 
chaque direction de traitement, 

" des moyens qui, pour le point courant, sont aptes a calculer Tecart 
inferieur des courbures des 4 connexes au point courant de la premiere 
image (1), 

- des moyens qui, pour le point courant, sont aptes a calculer, dans la 
deuxieme image (2), son point homologue interpole par la carte des 
ecarts de position, 

- des moyens qui, pour le point courant, sont aptes a calculer Tecart entre 
sa courbure et celle du point homologue, 

- des moyens qui, pour le point courant, sont aptes a rejeter le point 
homologue si cet ecart est plus important que Tecart inferieur calcule sur 
la premiere image (1). 

10. Precede d'acquisition et d'appariement d'un couple stereoscopique 
d'images, caracterise en ce qu'il comprend I'acquisition d'un couple 
stereoscopique damages (1, 2) par un instrument d'acquisition unique (30) 
comprenant deux capteurs CCD (31, 32) situes dans le plan focal optique 
(19*), le premier capteur (31) permettant Tacquisition d'une premiere image 
(1) du couple, et le deuxieme capteur (32) permettant I'acquisition d'une 
deuxieme image (2) du couple. 

11. Precede selon la revendication 10, caracterise en ce qu'il comporte en outre 
les phases consistant a : 
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- traiter gr§ce a des moyens de traitement les deux images (1, 2) du 
couple stereoscopique, selon au moins une direction de traitement et 
selon differents niveaux de resolution du plus grossier au plus fin, en, 
pour chaque changement de niveau de resolution ; 

- determinant, en chaque point de la premiere image (1), une 
fenetre de correlation (3) optimale, 

- calculant, par correlation avec les fenetres de correlation 
determinees, les ecarts de position entre chaque point (40) de la 
premiere image (1) et son homologue radlometrlque (41) dans la 
deuxieme image (2), 

- obtenant, pour le niveau de resolution courant, une carte des 
ecarts de position entre chaque point (40) de la premiere image et 
son homologue radiometrique (41) dans la deuxieme Image. 

- effectuant un recalage barycentrique des points de ladlte carte 
obtenue. 

- obtenir, pour le niveau de resolution le plus fin, la carte des hearts de 
position entre chaque point de la premiere image (1) et son homologue 
radiometrique dans la deuxieme image (2). 

12. Precede selon la revendication 11, caracterise en ce qu'il comporte une 
phase conslstant a traiter grace a des moyens de traitement les deux images 
(1, 2) du couple stereoscopique selon une direction epipolaire. 

13. Procede selon la revendication 11, caracterise en ce qu'il comporte une 
phase consistant a traiter grace a des moyens de traitement les deux images 
(1, 2) du couple stereoscopique selon une direction epipolaire et une 
direction orthogonale a la direction epipolaire. 
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14. Procede selon Tune quelconque des revendications 11 a 13, caracterise en 
ce qu'il comprend en outre une phase de saisie par un utiiisateur d'au moins 
un parametre d'exhaustivite correspondant au choix d*un rapport signal sur 
bruit de correlation suivant chaque direction de traitement. 

15. Precede selon I'une des revendications 11 a 14, caracterise en ce qu'il 
comporte en outre une phase d'acquisition de donnees comprenant 
Tacquisition : 

- d'une image de bruit par chaque capteur CCD, 

- du type de la fenetre de correlation, 

- d'une direction epipolaire pour chaque image (1, 2) du couple pour un 
coefficient stereoscopique moyen, 

- d'au moins une valeur de Pecart de position maximum suivant chaque 
direction de traitement, 

16. Precede selon la revendication 15, caracterise en ce que le precede 
comprend en outre une phase de pretraitement des donnees consistant a : 

- caiculer un nombre de nfveaux de resolution suivant chaque direction 
de traitement en fonction de chaque valeur de Tecart de position 
maximum, 

- filtrer les images (1, 2) du couple stereoscopique d'images par 
convolution avec une fenetre de convolution, 

- filtrer Timage de bruit par convolution avec ladite fenetre de convolution, 

- interpoler dans une geometrie epipolaire le couple stereoscopique 
d'images. 

17. Precede selon Tune des revendications 11 a 16, caracterise en ce que la 
phase de traitement des images comprend en outre les etapes consistant a : 
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- effectuer une convolution des images (1, 2) du couple par une fenetre 
de convolution, 

- effectuer un zoonn de tallle 2 des images (1, 2) du couple 
stereoscopique, 

18. Precede selon la revendication 14, caracterise en ce que Tetape de 
determination de la fenetre de correlation optimale consiste a, pour chaque 
point de la premiere image : 

- oalculer une courbure de la premiere image pour chaque taille de 
fenetre de correlation (3), 

- calculer une courbure de Umage de bruit pour chaque taille de fenetre 
de correlation (3), 

- calculer, pour chaque taille de fenetre de correlation, le rapport (SNRc) 
des courbures precedentes, 

et a choisir la plus petite taille de fenetre de correlation telle que ledit rapport soit 
superieur a chaque parametre d'exhaustivite saisi par I'utilisateur, 

19. Precede selon Tune des revendications 11 a 18, caracterise en ce que la 
phase de traitement des images comprend en outre une etape conslstant a 
effectuer une interpolation de la deuxleme image (2) grace a une carte des 
ecarts de position obtenue pour un niveau de resolution precedent. 

20. Precede selon la revendication 19, caracterise en ce que la phase de 
traitement des images comprend en outre les etapes consistant a : 

- rejeter les points homologues issus de fenetres de correlation plus 
grandes et contenant des fenetres de correlation plus petites. 

- rejeter fes points homologues (40, 41) de rayons de courbure verifiant 
une condition de rejet apres interpolation de la deuxieme image (2). 
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- completer iterativement, dans la carte des ecarts de position 
correspondant au niveau de resolution courant, des points d'ecarts de 
position lacunaires en fonction des points d'ecarts de position existants. 

- lisser la carte des ecarts de position correspondant au niveau de 
resolution courant par une fenetre de convolution. 

- composer la carte des ecarts de position calculee pour le niveau de 
resolution courant avec la carte des ecarts de position calculee pour le 
precedent niveau de resolution afin de mettre a jour ladite carte des 
ecarts de position calculee pour le niveau de resolution precedent. 

21. Procede selon Tune des revendications 20, caracterise en ce que Tetape de 
rejet des points homologues verlfiant une condition de rejet consiste a, pour 
le point courant : 

- calculer sur la premiere image (1) les courbures de ses deux points 
voisins suivant chaque direction de traitement, 

- calculer Tecart inferieur des courbures des 4 connexes au point courant 
de la premiere image, 

- calculer, dans la deuxieme image, son point homologue interpole par la 
carte des ecarts de position, 

- calculer I'ecart entre sa courbure et ceile du point homologue, 

- rejeter le point homologue si cet ecart est plus important que Tecart 
inferieur calcule sur la premiere image. 

22. Precede selon Tune des revendications 11 a 21, caracterise en ce que la 
phase d'obtention, pour !e niveau de resolution le plus fin, de la carte des 
ecarts de position, comprend Tobtention : 

- d'au moins un tableau representant la carte des ecarts de position 
suivant chaque direction de traitement, 
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- d'au moins un tableau representant une carte des tallies des fen§tres de 
correlation employees en tout point non rejete de la premidre Image, 
suivant chaque direction de traitement. 



23. Precede selon I'une des revendications 11 a 22, oaracterise en ce que les 
fenetres de convolution et de correlation utilisees sont des fenetres de 
convolution et de correlation de type prolate. 



recue le 12/05/04 





nil 




41 



4 ' c 



FIG-1 



regue le 12/05/04 




2 / 9 



8 



PiBliiliiiiiii 
liiiiiiilSir'iiifi 



Wmm 



Bif^^ir jar jfif m 




r 



V 



9 



10 



-12 

13 



FIG.2 



1 er depot 



Modifi^e le 12/05/04 




regue le 12/05/04 

3/9 




FIG.3 




FiG.6 



Modifi^e le 12/05/04 




regue le 12/05/04 



5/9 



k-- ^^^^^ -^'.-q^J 



21 



20 




22 








3' 



--7-1 



.7 




^ -m^' -^dv : 



23 




Modifieele 12/05/04 



.2,00 



1,50 



Q. 

s 

o 



1,00 



0.50 



0,00 



Variation de la precision en altitude en fonction du B/H 

1^ 



0,000 



3.1 



0,050 





















•'•-.5, • 


'••'■•'/ 




















-Vr*'---''- 


^-i 


















'..m-.r 








••'.•''.rJr.f 






































"CVn . * 
























































































































::;')':(>" 






















































































































r'tVf.'M-"" 








/,"<<••"•;■ 




















-.nr:.,'".:r 




n u;} 








































































































































i.- ri"t 


























































.■'<:^"'.c,• 








































































































Vf.«r ! 
















































V'Jjt'!;v 




































































































































SW'-;''?-.'. 




































^-■^'.«,'' 




















































■f'>',{'i 














f'.l.iM.V 




















































Vii'i:.': 


-UHi'Si 
































tAl 
















































































































































































































































































































































nViaHh 




















































































m 






~^ 






































it .th-'i 




















































































































































(v.. 














































































































































?.'!*?,^^,^ 


i-crjv 




































































hX-y-tii 






















'-•!>' A? ^ 
































































































































-4.r.'s-.<r 




































































































































































































fK-- 






































































t^iV 
















































































}iV-l*rl-l 
























































i.t.m'(-i' 

























































































































































0,100 



0,150 

B/H 



0,200 



0,250 



0,300 



Correlation Standard Prolate 
correlation fine 



—a— Correlation Standard Porte 

—o— Correlation fine Sans Inclusion de Fenetres 




6/9 



QJ 




rvT T-^ T-^ cT 



regue le 12/05/04 



7/9 




FI6.9 




F1G.10 



1er depot 



Modifies le 12/05/04 



) 

1. 




7. 



3oL 



^03 



:9 



M ml kst^ <^:^& 





1® 








M 1 






J .IJ 







f^: ^^^^ 




regue le 12/05/04 



9 / 9 



1 


RI1 




RI2 




RB 


1 

1 



301 302 303 

r i 



. 1 
1 


rh 




RI2 




RI3 


1 1 
1 i 



L. 



<j ? &» 



7 

304 



3 



r 



t 
1 


RH 




RI2 




RB 


1 
1 



1 
1 


RI1 




RI2 




RB 


T 

1 



300 



r 





RI1 




RI2 




RB 






1 




RI1 




RI2 


1 
1 


RB 





I 

10 



320 



r 



. 1 

1 


RI1 




RI2 




RB 


1 
1 





RI1 




RI2 


( 


R13 





330 



r 



1 


RH 




RI2 




RB 


1 

1 





RI1 


1 
1 


RI2 




RB 





J 



340 



FIGJ1 



KAYIDHAl OC 
L* FROPNICTe 
INDUtTRICLLr 



26 bis, rue de Saint Petersbourg - 75800 Paris Cedex 08 
Pour vous informer : INPI DIRECT 



0 82S 83 85 ST) 



0,15 € ITWimi 

Telecopie : 33 (0)1 53 04 52 65 



regue le 21/04/04 

lBi?E¥ET 0*!i^¥EP!lTI©M 
CERTIFICAT D'UTILITE n- 11235 03 

Code de la propriete intellectueMe - Livre VI 

DESIGi^ATION D'INVENTEUR(S) Page N° .1 . A - . . 

(A fournir dans le cas oD les demandeurs et les 
inventeurs ne sont pas les memes personnes) 

Get imprfmd est a remplir iisiblement a I'encre noire db 113 @w/ 210103 



my 



Vos references pour ce dossier (facultatij) 



TITRE DE L'INVENTION (200 caracteres ou espaces maximum) 

APPARIEMENT FIN D'IMAGES STEREOSCOPIQUES ET INSTRXJMENT DEDIE AVEC UN 
COEEFICIBNT STEREOSCOPIQTJE. 



FAIBLE 



LE(S) DEiyiANDEUR(S} : 

CENTRE NATIONAL D'ETUDES SPATIALES : 2, place Maurice Quentin 75001 PARIS FRANCE - 



DESIGI^E(NT} EN TANT QU'[nVENTEUR(S) 



in Nom 


ROUGE Bemai-d 


Prenoms 




Adresse 


Rue 


1% me Bonnat 31400 TOULOUSE FRANCE 


Code postal et vilfe 


1 t 1 ! 1 1 


Societe d'appartenance (facultatif) 




a Nom 




Prenoms 




Adresse 


Rue 




Code postal et ville 


1 1 I 1 1 1 


Societe d'appartenance (facultalif) 








Prenoms 




Adresse 


Rue 




Code postal et ville 


1 t 1 i i 1 


Societe d'appartenance (facultMf) 





S'ii y a plus de trois inventeurs, utflisez plusieurs formulaires. fndiquez en haut a droite le W de la page suivi du nombre de pages. 



DATE ET SIGNATURE(S) 
DU (DES) DERSANDEUR(S) 
OU DU EVIANDATAIRE 
(Nom et qualite du signataire) 




La ioi n^78-17 du 6 janvier 1978- relative a I'lnformatique, aux ftchiers et aux libertes s'applique aux reponses faites a ce formulaire. 
Elle garantit un droit d'acces et de rectification pour les donnees vous concernant aupres de TINPL 



